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Abstract. The interference between coherent and squeezed vacuum hght can 
produce path entangled states with very high fidelities. We show that the phase 
sensitivity of the above interferometric scheme with parity detection saturates the 
quantum Cramer-Rao bound, which reaches the Heisenberg-limit when the coherent 
and squeezed vacuum light are mixed in roughly equal proportions. For the same 
interferometric scheme, we draw a detailed comparison between parity detection and 
a symmetric-logarithmic-derivative-based detection scheme suggested by Ono and 
Hofmann. 
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1. Introduction 

Optical metrology relies on light interferometry as its primary tool for phase estimation. 
The sensitivity of phase estimation with coherent light based interferometry is limited 
by shot noise [1]. This limit, however, is due to the classical nature of coherent light and 
can be surpassed if nonclassical states of light, such as the NOON state, are used [21 [3]. 
Still there is a limit on the sensitivity of phase estimation in the case of linear optical 
interferometry. Its usual justification stems from the Heisenberg uncertainty principle 
that links phase uncertainty of a state to its photon number uncertainty, A(j)An > 1. 
Combination of this equation with the assumption that the photon number uncertainty 
in a state is limited by the total photon number (in the case of states with definite photon 
number) or the total average photon number (in the case of states with indefinite photon 
number). An < A^, suggests the limiting phase sensitivity to be A(/)hl = which is 
commonly referred to as the Heisenberg limit [H E] . 

Quantum optical metrology has Heisenberg limited sensitivity of phase estimation 
as its goal. To this end, search for convenient states of light and optimal detection 
schemes still continues P, El [HI [9l [TOl [11]. Candidate states of light are gauged based 
on the quantum Cramer-Rao bound that provides a detection scheme independent 
phase sensitivity A(/)qcrb [E]. In turn, optimal detection schemes are sought, which 
are capable of saturating the quantum Cramer-Rao bound. A known possibility is a 
detection scheme that measures a symmetric logarithmic derivative, since such operators 
saturate the quantum Cramer-Rao bound; however, they are seldom easy to implement. 
The capabilities of alternative detection schemes are judged by the classical Cramer- Rao 
bound that is detection scheme specific, or by the error propagation formula that links 
the uncertainty of the observed signal with phase uncertainty. 

Here, we consider coherent and squeezed vacuum light input as a candidate for 
Heisenberg limited phase estimation (see Fig. [1]). This state has been previously 
checked against the quantum Cramer-Rao bound and shown to achieve Heisenberg 
limited phase sensitivity when the coherent and squeezed vacuum light are mixed in 
roughly equal intensities — a feature that can be explained as due to the high fraction 
of a NOON state in the normalized N-photon output component of the quantum state 
of light past the mixing beam splitter — first pointed out by Hofmann and Ono [14j and 
later experimentally demonstrated by Silberberg's group [I5]. The detection scheme 
suggested in Ref. [13j for Heisenberg limited phase estimation was however based on 
Bayesian analysis of the photon number statistics of the output state, which requires 
photon number counting in both modes. 

In this paper, we study parity detection p£j for the interferometry with coherent 
and squeezed vacuum light. We show that the parity operator saturates the quantum 
Cramer-Rao bound, and in turn provides Heisenberg limited phase sensitivity when 
the coherent and squeezed vacuum light are mixed in equal proportions. Parity 
detection should be a simpler alternative to the detection scheme of Ref. [13] since 
parity measurement can be inferred from the photon number counting statistics of a 
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Figure 1: A lossless Mach-Zehnder Interferometer with two-mode state input given by 
the product of coherent and squeezed vacuum states of light. The modes are labeled 
by the annihilation operators: a^, hi for the input ports and a/, bf for the final output 
ports, respectively. BSl and BS2 are 50-50 beam splitters. 

single mode alone, which in the low power regime can be obtained using phot on- number- 
resolving detectors [17J. Although accurate photon-number-resolution in the high power 
regime is a difficult task, there have been proposals for the quantum non-demolition 
measurement of photon number using weak non-linearities and homodyning [18j. On 
the other hand, it is not necessary to have phot on- number-resolving capabilities in order 
to implement parity detection. Assuming the availability of large Kerr nonlinearities 
through the techniques of electromagnetically induced transparency, a scheme that 
performs quantum non-demolition measurement of parity directly, without requiring 
the measurement of photon number, has been proposed [I9]. Plick et al have recently 
shown that parity measurement for interferometric schemes that use Gaussian states, 
like the one in use here, could possibly be inferred through balanced homodyning and 
intensity difference measurement [20] . 

Ono and Hofmann, in Ref. [21], have studied a detection scheme based on the 
measurement of a symmetric logarithmic derivative for the considered interferometric 
scheme. This scheme uses interference with an auxiliary local oscillator and intensity 
difference measurement as well. We dutifully discuss this scheme with an intent to 
compare it with parity detection. 

The paper is organized as follows. Section [2] describes the propagation of a two mode 
light, initially in the product state of coherent and squeezed vacuum light, through the 
Mach-Zehnder interferometer. Section [3] focuses on the parity-based detection scheme 
and provides the expected signal and phase sensitivity, while section [H discusses the 
Ono-Hofmann detection scheme in equal detail. Section [5] deals with the conclusion. 

2. Propagation of the input fields through the Interferometer 

The input to the interferometer is in the product state |c^o) ® |C = ^ ^^'^^) that describes 
coherent light with displacement ao = y^e"^"^^ in one mode and squeezed vacuum with 
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parameters r and 0^ in the other. The corresponding Wigner function of the input state 
is the product of the respective Wigner functions as weU [22] : 

WUcy, ao; /3, r) = W,{a, ao)Ws{/3, r), (1) 

with Wigner function for the corresponding states being 

W,{a,ao) = fe-2|"-"oP, WsiP,r) = le-2|/3pcosh2.-(/32+/3-)sinh2.^ (2) 

and where we have made = by appropriately fixing the irrelevant absolute phase. 
This choice implies that the phase of the coherent light (/)c is now measured with respect 
to the phase of the squeezed vacuum state. 

A Mach-Zehnder interferometer is composed of optical elements such as beam 
splitters, mirrors and phase shifters. Propagation of the light field through these 
elements is described by relating the initial variables in the Wigner function to their 
final expressions: 

Wouticxf, f3f) = WUc^,{af, I3f),/3,{af, /?;)). (3) 

The relation between variables in the most general form is given by a two-by-two 
scattering matrix M: 

(4) 

where a^, a/, and /3/ represent the complex amplitudes of the field in the modes a^, 
bi^ a/, and 6/, respectively. More specifically, propagation through a 50-50 beam-splitter 
and a phase shifter (in mode 6), are described by: 
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respectively. Therefore, the Mach-Zehnder interferometer in Fig. [T] is described by 
Mmzi = MbsM^Mbs and is found to be: 
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with the corresponding transformation of the variables in the following form: 

ai —le 2 [^OLf sm — + cos — j, 

A -ie'^af cos^ - f3f sin^). 

Therefore, the state of light at the output of the Mach-Zehnder interferometer is 
described by the following Wigner function: 



(6) 



(7) 
(8) 
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Having found the state of light at the output of the Mach-Zehnder interferometer, 
we will present the parity-based phase estimation scheme with calculations of its signal 
and phase sensitivity in the following section. 
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3. Phase Estimation with Parity Measurement 

Parity detection was originally proposed in the context of trapped ions by Bollinger et 
al [23]. It was later adopted for optical interferometry by Gerry [24] . In its essence, 
parity detection distinguishes states with odd and even number of photons. Quantum 
mechanically it is described by the parity operator, = (—1)^^^, acting on a single 
output mode, a. Parity detection makes phase inference at the Heisenberg limit possible 
without having to know the full photon number counting statistics for several classes of 
input states with definite as well as indefinite photon numbers (including the NOON 
state) [25l [26]. Coherent and squeezed vacuum light belong to the latter class of 
states and the performance of parity detection for these states is studied in this section. 
Although the fact that parity detection achieves Heisenberg limited performance with 
the NOON-state is a motivation for this study, it is by no means a reason in itself to believe 
that parity detection would work equally well for the considered interferometric scheme. 
This is because the quantum state in the interferometer also has several miscellaneous 
contributions apart from that of a NOON state. 

An expected signal of the parity detection scheme (Ha) is calculated as the value 
of the Wigner function at the origin for the corresponding mode. In the case of mode 
a/, (Hay) = 2 / ^out(0, and is found to be: 



exp 
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(9) 
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where the coherent light amplitude and the squeezing parameter have been expressed 
in terms of the average photon numbers, Uc and n^, using the relations ao = y^e"^"^^ 
and r = sinh""^ 

The signal of the parity detection scheme is periodic with period 27r and attains its 
maximum value of one at = 0. Although this maximum value is independent of the 
phase of the coherent light (j)c and the light intensities Uc and n^, the visibility of the 
signal and its width are functions of these parameters. The visibility of the signal is found 
to be best when (j)c = and to diminish as 0c drifts away from zero, becoming worst 
at 0c = 7r/2. Since it is reasonable to assume the coherent and squeezed vacuum light 
to be locked to the same external phase, 0c can be set to zero for optimal performance. 
Further, the dependence of the signal on the light intensities is studied in terms of the 
total input intensity, riin = nc + n^, and the fraction of total intensity in the squeezed 
vacuum state, rj = Ug/nin- When rj is increased from zero, the signal is found to grow 
narrower until reaching an optimal width, and then to broaden again, but with reduced 
visibility as r] approaches one. For rj = and = 1, the width of the signal is found 
to be proportional to Tr/y^ni^, which is narrower than the resolution of conventional 
interferometry by a factor of y^r?^ and thus demonstrates super-resolution [27] . The 
fraction 77 = 0.5 is found to be the most optimal choice for distributing the input 
light intensity, since it allows a higher narrowing factor of Ui^. Figure [2] demonstrates 
this result by comparing the parity signals for interferometry with only coherent light 
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(77 = 0) [28] or squeezed light (ry = 1) and interferometry with coherent and squeezed 
vacuum hght of equal intensities (ry = 0.5). We see that for the same total input photon 
10, the parity signal for the latter case is narrower than any other case. 



number, 




-7r/2 -7r/3 -7r/6 



7r/6 7r/3 7r/2 



Figure 2: (Color online) The parity signal, (Ha/), as a function of the accumulated 
phase difference between the arms of the MZI, 0: dashed (magenta) line for coherent 
light interferometry (77 = 0) with Uc = 10; dotted (red) line for squeezed vacuum light 
interferometry (ry = 1) with Ug = 10; and solid (blue) line for coherent and squeezed 
vacuum light interferometry {rj = 0.5) with Uc = Ug = 5, (pc = 0. The dot-dashed 
(green) line is the signal for conventional coherent light interferometry with intensity 
difference measurement. 



The phase sensitivity A(j) of an interferometer, followed by a detection scheme 
described by an operator O, can be characterized using the error propagation formula: 

id') - {Of 



d{d)/d4> 



(10) 



For the parity based detection scheme, O 



^af> knowing the signal suffices for 
sensitivity calculation since H^^ = 1. The phase sensitivity with parity detection for 
coherent and squeezed vacuum light interferometry is found to be best at = 0, and is 
given by: 

A^' = z . . ^— — • (11) 



2nc\/ ns{ns + l)cos 2(/)c + IncUs + nc + ris 
For a detection scheme to be optimal, it has to saturate the quantum Cramer-Rao 
bound. The quantum Cramer-Rao bound for the considered interferometric scheme was 
derived in Ref. [13] and found to be: 

1 

^QCRB 
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(12) 



This expression can be shown to be identical to the phase sensitivity with parity 
detection, given in Eq. ( ITTI) (under the condition (/)c = 0), when and r are replaced in 
terms of the average photon numbers, Uc and rig. Thus, parity detection saturates the 
quantum Cramer-Rao bound and is optimal for the considered interferometric scheme. 
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Figure 3: (Color online) Phase sensitivity with parity detection for coherent and 
squeezed vacuum light interferometry A(/) as a function of the fraction of squeezed 
vacuum in the input r]. The total input photon number = 10 and (/)c = 0. 
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Figure 4: (Color online) Phase sensitivity with parity detection A(/) as a function of the 
accumulated phase difference between the arms of the MZI dashed (magenta) line 
for coherent light interferometry (77 = 0) with Uc = 10, (j)c = 0, dotted (red) line for 
squeezed vacuum interferometry (ry = 1) with = 10 and solid (blue) line for coherent 
and squeezed vacuum light interferometry (77 = 0.5) with nc = = 5, (/)c = 0. 

Although parity detection is optimal for the considered interferometric scheme 
irrespective of the input intensities, the combination as a whole achieves its best phase 
sensitivity when = 0.5. Figure [3] is a plot of the phase sensitivity Acj) given in Eq. f fTTI) 
(under the condition (/)c = 0), as a function of the fraction of squeezed vacuum in the 
input T]. The phase sensitivity can be seen to be best when 0.5. Equation f fTTI) . under 
the condition (/)c = 0, reveals that the phase sensitivity of the combination coincides with 
the Heisenberg-limit, Acj) ^ l/rzin, when 77 ?^ 0.5, while it coincides with the shot-noise 
limit, A(j) ^ l/y^rii^, when = or 1. Figure [4] illustrates the contrast in phase 
sensitivities with parity detection for the cases corresponding to coherent and squeezed 
vacuum light interferometry with = 0, 77 = 1, and = 0.5. 

Finally, we study the critical dependence of phase sensitivity on the fraction of 
squeezed vacuum r] and the input phase 0c near the optimal point of operation, namely 
T] = 0.5, (j)c = 0. The contour plot given in Fig. [5] shows that with small variations in 77 
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Figure 5: (Color online) Contour plot of the phase sensitivity with parity detection for 
coherent and squeezed vacuum light interferometry close to phase origin, as a function 
of the phase of the coherent light (^c and the fraction of squeezed vacuum in the input 
77, respectively. The total photon number = 10. At the point 77 = 0.5 and 0c = 0, 
phase sensitivity attains it's optimum value of 0.093. 

and (/)ci the phase sensitivity of the scheme is not much compromised. Thus, the parity 
detection scheme provides Heisenberg limited phase sensitivity even while allowing for 
small fluctuations in rj and (j)c about the optimal point of operation. 

4. Discussion 

So far, we have shown that parity detection could be used to achieve Heisenberg- 
limited phase estimation in the interferometry with coherent and squeezed vacuum 
light. The different possible implementations of parity detection include the use of 
phot on- number-resolving detectors in the low power regime; Kerr nonlinearities and 
homodyning in the high power regime. In Ref. [21j, Ono and Hofmann discussed a 
different detection scheme that implements the measurement of a symmetric logarithmic 
derivative. Implementation of this measurement is based on interference with a local 
oscillator and intensity difference measurement (see Fig. [6]). Since symmetric logarithmic 
derivative based phase estimators saturate the quantum Cramer- Rao bound, Heisenberg- 
limited phase sensitivity was anticipated with this scheme for the interferometry with 
coherent and squeezed vacuum light mixed in equal proportions (77 = 0.5). Here, we 
present a brief study of the Ono-Hofmann detection scheme (in the absence of losses), 
for the purpose of comparing it with parity detection. 

The Ono-Hofmann detection scheme consists of a second MZI appended at the 
output of the flrst, with a control phase 99, which is set to zero[j|. A local oscillator fleld, 

J It was set to tt by Ono and Hofmann in their original analysis. We choose zero for ease of calculation. 




[>QlBl 
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which is in the coherent state, I710) = A/nio/Te^^^^^, is introduced by mixing with the 
mode CLff through a highly reflective beam sphtter of transmissivity, T << 1, where nio 
is the average number of photons in the fleld that eventually enters the interferometer, 
and 010, its phase. In the end, the difference in intensities at the two output modes is 
measured. 
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Figure 6: The Ono-Hofmann detection scheme for interferometry with coherent and 
squeezed vacuum light. The detection scheme uses interference with an auxiliary local 
oscillator and intensity difference measurement for phase estimation. A highly reflective 
beam splitter is used to mix the local oscillator fleld into the interferometer. 



Intensity measurements at the output provide: 

(13) 



(CO = ({C^c}^/ 2' 



{c^c}^ ) being the symmetric form of the operator, which can be evaluated 

based on the flnal Wigner function of the state Wf as: 

({c^c},) = J J \a\^Wf{a,l3)d^ad^l3, (14) 
(jrftrf} ) = J J \P\^Wf{a,p)d^ad^P, (15) 

where a and /3 are the complex amplitudes in the modes c and d respectively. 

The signal, which is the difference in intensities at the output ports, is thus given 

by: 

/ = / / (lap - |/3|2) Wf{a, /3)d^ad^/3, (16) 
and is found to be: 

/ = —2^ncn\o cos | cos (f + 0c — 0io) + (^c — ^s) sin^. (17) 

It is plotted in Fig. [71 as a function of 0, under the condition 0c = 0, 0io = 7r/2 (the point 
where the phase sensitivity is found to be optimal, as mentioned later in the paper). 

The phase sensitivities in both cases turn out to be identical, although the signals differ. 
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The figure compares the signal for the interferometry with equal intensities of coherent 
and squeezed vacuum light (77 = 0.5), with those of interferometry with only coherent 
light {t] = 0) and only squeezed vacuum light (77 = 1). We see that for the same total 
input photon number, = 10, and same number of photons from the local oscillator, 
nio = 100, the signal for the former is stronger than any other case. However, unlike 
with parity detection, there is no super-resolution in the signal for the Ono-Hofmann 
detection scheme. 
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Figure 7: (Color online) The signal with the Ono-Hofmann detection scheme — the 
intensity difference / is plotted as a function of the accumulated phase between the 
arms of the MZI 0: dashed (magenta) line for coherent light interferometry (77 = 0) with 
TT^c = 10, 0c = 0, dotted (red) line for squeezed vacuum light interferometry (77 = 1) 
with 7^5 = 10 and solid (blue) line for coherent and squeezed vacuum light interferometry 
(77 = 0.5) with He = Us = 5, (j)c = 0. A local oscillator of strength 77.10 = 100 and phase 
0io = is used. 




Figure 8: (Color online) Phase sensitivity of coherent and squeezed vacuum light 
interferometry as a function of the fraction of squeezed vacuum light in the input ij 
^in = 10: solid (magenta) line for phase sensitivity with the Ono-Hofmann detection 
scheme, A(/), dashed (blue) for the Quantum Cramer- Rao bound, /^(^qcrb- 

Phase sensitivity with the Ono-Hofmann detection scheme for the interferometry 
with coherent and squeezed vacuum light is calculated based on the error propagation 
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formula mentioned in Eq. (fTOl) . Variance of the signal AP which is required in the 
formula, can be shown to be: 



1 

2' 



The phase sensitivity thus calculated, is found to be optimal at 



condition 



0, 



^lo = 7r/2, and given by: 

. _ -2(-2ns+2v/^^~fTv^-l)(v^-V^)2+2ns+l 



(18) 
TT, under the 



(19) 



In the limit of infinite nio, it takes the following simplified form: 



^02 2ns-2^yn^^+ns+l ^ 1 



for {us > 1). 



(20) 




Figure 9: (Color online) The Ono-Hofmann detection scheme for the interferometry with 
coherent and squeezed vacuum light: A contour plot of A(j) x riin with respect to and 
the fraction of squeezed vacuum light in the input r] in the presence of a local oscillator 
of infinite power. The contour of value one corresponds to points of Heisenberg-limited 
phase sensitivity. It can be seen that this contour tends to 77 = 0.5 in the limit of large 



This expression nearly saturates the quantum Cramer-Rao bound for the 
interferometry with coherent and squeezed vacuum light given in Eq. f fT2l ). as can be 
seen in Fig. [HI which compares graphs of the quantum Cramer Rao bound and the 
phase sensitivity given in Eq. f l20l ). as a function of the fraction of total intensity in the 
squeezed vacuum state rj for a given total input intensity, = 10. The fact that they 
don't completely overlap can be explained as due to the use of the error propagation 
formula instead of the more accurate classical Cramer Rao bound, which in this case is 
difficult to calculate. 

Since we are mainly interested in achieving Heisenberg limited phase sensitivity 
with the considered interferometric scheme, we now focus our attention on the specific 
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Figure 10: (Color online) The Ono-Hofmann detection scheme for the interferometry 
with coherent and squeezed vacuum light: A contour plot of A(j) with respect to 

and the fraction of squeezed vacuum light in the input rj in the presence of a local 
oscillator with photon number nio = 100. A contour of value one would correspond to 
points of shot noise-limited phase sensitivity, while contours of values above one refer 
to points where phase sensitivity is worse than shot noise. 



case corresponding to 77 = 0.5, since the quantum Cramer Rao bound of Eq. f fT2l ) for 
this 77 matches the Heisenberg limit. Having found that the Ono-Hofmann detection 
scheme saturates the quantum Cramer Rao bound (in the limit of infinite nio), we look 
to estimate the amount of local oscillator power that is required to provide Heisenberg- 
limited phase estimation, for a given total input photon number, n. In other words, 
we look to make an order of magnitude estimate for the infinity in terms of total input 
intensity. When Ug ^ Uc ^ the expression in Eq. (fT9l) takes the form: 

^2 / (^in + l- V^hi + 2^in) (nin + 2nio-V8^inrilo)+^in + l J^2l) 

In the limit of large nin, the regime of interest of the Ono-Hofmann detection scheme, 
this can be expanded in a series as: 

1 S 1 

= ^= + . V- - ^= (22) 

From the above expansion, it is evident that the local oscillator photon number, nio, has 
to be as large as n?^ for the scheme to provide Heisenberg-limited phase sensitivity. Thus, 
if one were to include the local oscillator in the photon budget, the phase sensitivity 
that the Ono-Hofmann detection scheme can be expected to provide is at best shot-noise 
limited. This shortcoming is illustrated using Figs. [9]and[ini which show contour plots 
drawn with respect to the total input photon number and the fraction of total input 
intensity in the squeezed vacuum state 77. Although the total photon number in the 
plots run only up to ten, they still illustrate the point. Figure [9] shows the plot of the 
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product of phase sensitivity of the Ono-Hofmann detection scheme in the presence of a 
local oscillator of infinite power, as given in Eq. f l20l) . and the total input photon number 
riin- The presence of contours of value close to one indicates that phase sensitivity with 
the Ono-Hofmann detection scheme is indeed Heisenberg-limited when the coherent and 
squeezed vacuum photons alone are counted as photon resource. Figure [IHl on the other 
hand, shows the product of phase sensitivity of the Ono-Hofmann scheme in the presence 
of a local oscillator of finite power, as given in Eq. ffTOl)). and the square root of the sum 
of photon numbers in the input and the local oscillator field, = ^/rv^F^^- It is, 
however, seen that there exist no contours of value one or less, which indicates that the 
phase sensitivity is in fact worse than shot noise when the photons in the local oscillator 
are also accounted for as part of the photon resource. 

At this point, we'd like to mention that the scheme proposed by Plick et al. in 
Ref. [20] to implement parity detection for interferometry with Gaussian states using 
balanced homodyning and intensity diflFerence measurement would also suflFer from 
the same drawback when the photons in the local oscillator are counted as resource. 
Nevertheless, given the fact that these schemes are proposed for applications in the 
high power regime, the above mentioned drawback is not so consequential. At LIGO, 
while there is plenty of laser power available, the interferometers are power-limited 
with the risk of melting the mirrors and beam-splitters. Therefore, the power in the 
local oscillator in any case cannot be used inside the interferometer. Thus, phase 
estimation, which is Heisenberg-limited in the photon number inside the interferometer, 
as provided by the Ono-Hofmann detection scheme or parity measurement inferred 
through homodyning, would still be appreciated. 

5. Summary 

We have studied the application of parity detection for phase estimation in Mach- 
Zehnder interferometry with coherent and squeezed vacuum light. We have shown 
that parity detection saturates the Quantum Cramer Rao bound of the interferometric 
scheme and provides Heisenberg-limited phase sensitivity when the coherent and 
squeezed vacuum light are mixed in equal proportions. Parity can be readily 
implemented using photon- number-resolving detectors [17j in the low power regime, and 
possibly using optical nonlinearities and homodying in the high power regime p!8l[T9l[20] . 

We have also presented a brief study of a symmetric-logarithmic-derivative-based 
detection scheme recently proposed by Ono and Hofmann for the same interferometric 
scheme in the high power regime [2T] , with explicit calculations of the signal and phase 
sensitivity of the scheme. We have shown that this scheme requires a strong local 
oscillator field in order to provide Heisenberg-limited phase sensitivity. When the local 
oscillator power is accounted for as part of the photon resource, the phase sensitivity 
of the scheme is at best short-noise limited — a drawback parity detection implemented 
using homodyning is also bound to suffer. 
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